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' Abstract. We present results from a systematic study of the effects of initial 

' parameters on three-dimensional thermonuclear supernova models. 

Oh. 

' One major challenge to numerical models of thermonuclear supernova ex- 

plosions is to reproduce the diversity observed among Type la supernovae (SNe 
la). With the rapid advancement of three-dimensional models over the past 
years, systematic studies have become possible testing the effect of changes in 
the initial parameters of simulations on the explosion process. The ultimate 
goal is, of course, to explain the correlation between the peak luminosity and 
I the light curve shape used to calibrate cosmological distance measurements. 

Tij" ■ We present the first systematic study on this issue based on three-dimen- 

0^ I sional deflagration models. Here, a Chandrasekhar-mass white dwarf (WD) star 

ignites nuclear reactions in the center which finally form a flame. In the defla- 
I gration model this flame burns outward with a subsonic velocity, which quickly 

becomes dominated by the effects of interaction with a turbulent velocity field. 
Caused by buoyancy (Rayleigh- Taylor) and secondary shear (Kelvin-Helmholtz) 
instabilities inherent in the scenario, turbulence wrinkles the surface of the flame 
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2 ! effectively boostin g its propagation speed. For a c omprehensive review of SN la 

' models we refer to iHillebrandt Sz Niemeveil (2003)- 

. For the three-dimensional simulations of the explosion process we apply the 

scheme developed by iReinecke et al. (1999) and Rcinecke et al. ( 2002i)- The final 
composition of the ejecta in our models is obtained via a nucleosynthetic post- 
^ ' processing procedure on the basis of tracer particles advected in the explosion 

■ simulation (for details see iTravag lio et al. , 20041. 

There is a number of possible parameters of the deflagration model that 
have the potential to explain the SN la diversity. The progenitor's carbon- 
to-oxygen (C/0) ratio, its metallicity, and the central density at ignition are 
commonly suggested, but other effects like rotation and flame ignition could 
also play a role. In our survey we concentrate on the former three parameters. 
These are varied independently. Of course, this is an oversimplification since 
in reality they are interrelated by stellar evolution of the progenitor WD star. 
Nevertheless, in this first study our goal is to infer the trends of effects of each 
individual parameter on the explosion. 

In the setup of the explosion models we apply three different carbon mass 
fractions of the WD material, X{^'^C) = 0.30,0.46,0.62, and three different 
central densities at ignition, pc = [1.0,2.6,4.2] x lO^gcm"^. Three different 
metallicities of the WD, Z = [0.3, 1.0, 3.0]^©, are represented by the ^^Ne mass 
fractions fsee lTimmes et al]l2003t ). This defines the 27 models of our survey. 
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A change in the carbon mass fraction does not alter the flame evolution in 
our models significantly. This is due to a larger amount of a-particles in the 
nuclear statistical eq uilibrium in the ashes for carbon-rich fuel, which act as 
an energy buffer fsee lRooke Sz Hillebrand"3l2004^ ■ This effect delays the flame 
propagation at explosion stages where iron group nuclei are synthesized and 
therefore the changes in the resulting ^^Ni masses are little. Since the radioactive 
decay of this isotope powers the lightcurve, according to Arnett's rule the peak 
luminosity is nearly unaffected by the C/0 ration of the progenitor. However, in 
the later evolution the carbon-rich fuel models release more energy (amounting 
to a ^12% variation) which possibly affects the light curve shape. 

A lower central density leads to a lower energy release and delays the evo- 
lution of the model. The reason for this effect is the reduced gravitational 
acceleration experienced the flame front which leads to a slower formation of 
the nonlinear Ray leigh- Taylor instabilities. We find a ~40% variation in the 
energy releases of our models. The amount of produced ^^Ni is largest for the 
intermediate central density. At lower densities the delayed flame propagation 
leads to a decreased production of iron group elements. In models with higher 
central densities electron captures favor neutron rich nuclei instead of ^^Ni. This 
effect is accounted for in the nucleosynthetic postprocessing but not consistently 
modeled in the explosion simulation. Thus the results are preliminary for the 
highest central density value. The change in the produced ^^Ni is ~10%. A 
variation in the central density prior to ignition is likely to affect both the peak 
luminosity and the light curve shape. 

Changing the metallicity (i.e. the ^^Ne mass fraction) we find a ~20% 
change in the resulting ^^Ni masses and h ence the peak lumino sitv of the event 
consistent with the analytic prediction of iTimrnes et al.l (|2003). The explosion 
dynamics, however, is not affected by the metallicity and thus the light curve 
shape will not differ significantly. 

Concluding we note that the variations in the peak luminosities obtained by 
changing the initial parameters can partially account for of the observed scatter 
in SNe la. However, we could not identify a single parameter that reproduces 
the empirically established peak luminosity-light curve shape relation. A combi- 
nation of initial parameters chosen according to stellar evolution could possibly 
explain such a relation, but initial parameters ignored in this first survey have 
to be explored as well in future studies. The final decision on the effect of the 
initial parameters on the light curve shape can only be made on the basis of 
synthetic light curve calculations. 
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